Multiwavelength study of the starburst galaxy NGC 7714. I: 
Ultraviolet- Optical spectroscopy 1 
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ABSTRACT 

We have studied the physical conditions in the central 300 pc of the proto-typical 
starburst galaxy NGC 7714. Our analysis is based on ultraviolet spectroscopy with the 
HST+GHRS and ground-based optical observations, and it covers also X-ray and radio 
data taken from the literature. The data are interpreted using evolutionary models 
optimized for young starburst regions. The massive stellar population is derived in a 
self-consistent way using the continuum and stellar absorption lines in the ultraviolet 
and the nebular emission line optical spectrum. 

The central starburst has an age of about 4.5 Myr, with little evidence for an age 
spread. Wolf-Rayet features at the ultraviolet indicates a stellar population of ~ 2000 
Wolf-Rayet stars. The overall properties of the newly formed stars are quite similar 
to those derived, e.g., in 30 Doradus. A standard Salpeter IMF is consistent with all 
observational constraints. The nucleus of NGC 7714 has a bolometric luminosity of 0.5- 
1 x 10 10 L© and a mass of 5-10 x 10 6 M© (if low-mass stars have formed) . We find evidence 
for spatial structure within the central 300 pc sampled. Therefore it is unlikely that the 
nucleus of NGC 7714 hosts a single star cluster exceeding the properties of other known 
clusters. Contrary to previous suggestions, we find no evidence for a nuclear supernova 
rate that would significantly exceed the total disk-integrated rate. About one supernova 
event per century is predicted. Most of them are associated with the core-collapse of a 
hydrogen-free or -poor progenitor. An older stellar generation, with ages of tens of Myr 
and older, is suggested as well. This population is less concentrated towards the nucleus 
and extends over kpc scales 

Subject headings: galaxies:evolution-galaxies:starburst-galaxies:stellar content-galaxies: ultraviolet- 
galaxies:ISM-stars:early-type-galaxies:individual:NGC 7714. 
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1. Introduction: NGC 7714 a test case for 
evolutionary synthesis models 

Starburst galaxies are objects in which the total 
energetics are dominated by star formation and asso- 
ciated phenomena (Weedman 1983; Heckman 1998). 
This definition covers galaxies with a very wide range 
of properties, from blue compact dwarfs to ultralu- 
minous IRAS starbursts. Typical masses (bolomct- 
ric luminosities) range from 10 6 to 10 10 M© (10 9 
to 10 14 Lo), corresponding the lowest limit to the 
mass of super-star-clusters and the highest limit to 
the mass of the infrared- luminous galaxies (Leitherer 
1996). Nearby starburst galaxies are ideal laborato- 
ries in which to explore fundamental questions about 
the local and global processes and effects of star for- 
mation. For example, is the initial mass function the 
same in starburst galaxies as in quiescent star forma- 
tion in the Milky Way? Is there a connection between 
starburst activity and active galactic nuclei? What 
can starbursts teach us about the process by which 
the first galaxies formed? 

Starburst are powered by massive stars. Massive 
stars (M> 10 M©) emit photons with energies of tens 
of eV which are absorbed and re-emitted in their stel- 
lar winds, producing ultraviolet resonance transitions. 
However, the stellar wind is optically thin to most 
of the ultraviolet photons, that can travel tens of pc 
from the star before they are absorbed and photoion- 
ize the surrounding interstellar medium. Then, this 
ionized gas cools down via an emission line spectrum. 
This is, in essence, the spectral dichotomy picture of 
a starburst galaxy: a nebular emission line spectrum 
at optical wavelengths and an absorption-line spec- 
trum at wavelengths shorter than the Balmer jump 
(Leitherer 1997). However, the presence of a star- 
burst affects the entire spectrum and not only the 
UV-optical spectrum. Thus, supernova remnants fol- 
lowing the explosion of massive stars generate X-ray 
and radio continuum emission; red supergiant stars 
emit strongly in the near-IR continuum; the mid-IR 
is dominated by line and continuum emission from 
warm dust within the star forming regions; the far- 
IR is thermal emission from dust heated by absorbing 
UV and optical light; electron accelerated by super- 
nova remnants create radio synchroton emission. A 
truly complete analysis of a starburst would account 
for all these features. 

Evolutionary synthesis is a powerful tool for under- 
standing starburst galaxies. This technique makes a 



prediction for the spectrum of a stellar population 
taking as a free parameter the star formation his- 
tory of the starburst (age, initial mass function, star 
formation rate, etc). If there is agreement between 
the observations and the predicted physical proper- 
ties (spectral energy distribution, colors, line profiles, 
etc), constraints on the star formation history (in 
timescale between 1 Myr to several Gyr) can be deter- 
mined (see e.g. Mas-Hesse & Kunth, 1991; Bruzual 
& Chariot, 1993; Leitherer & Heckman, 1995). How- 
ever, the solution obtained with this technique may 
not be unique, and only consistency between models 
and observations can be derived. 

To make a complete analysis of a starburst, we 
have to answer the question: are evolutionary syn- 
thesis models, and the stellar evolution models they 
incorporate, actually up to the challenge of interpret- 
ing multi-wavelength data? To answer this, we have 
undertaken a detailed study of the starburst galaxy 
NGC 7714. 

NGC 7714 is an SBb peculiar galaxy (dc Vau- 
couleurs et al 1991) classified by Weedman et al (1981; 
here after W81) as a prototypical starburst galaxy. 
Its heliocentric radial velocity, as compiled by NASA 
Extragalactic Database (NED), is 2798 km s _1 . This 
places NGC 7714 at a distance of 37.3 Mpc assuming 
H = 75 km s^ 1 Mpc" 1 . Heckman et al (1998) de- 
rive a distance of 37.9 Mpc using the standard linear 
Virgo centric infall model of Schechter (1980) with 
parameters 7=2, v\/irgo= 976 km s _1 , z«0=22O km 
s _1 and Dvirgo= 15-9 Mpc. The B magnitud and 
far-infrared luminosity (as calculated from the 60 /im 
and 100 ^m IRAS bands) are -20.04 and 43.93 erg s" 1 , 
respectively. The far-infrared ratio 60/100 and 25/60 
are larger than 0.9 and 0.25, respectively. The mass 
of the neutral gas is ~ 10 10 M© (Mirabel & Wilson, 
1989). The X-ray, IUE ultraviolet, optical, and VLA 
radio fluxes are explained as a consequence of intense 
star formation activity in its nucleus. To explain the 
X-ray and radio luminosities, 10 supernova remnants 
are required (W81). Bcrnlohr (1993) found from the 
optical continuum and FIR colors that the star for- 
mation in NGC 7714 is consistent with a continu- 
ous star formation rate during the past 20 Myr, and 
that it had probably been triggered by the interaction 
with the companion galaxy NGC 7715 (Vorontsov- 
Velyaminov 1977; Arp 1966; Smith, Struck, & Pogge, 
1997). Very detailed studies have been carried out in 
the optical and near-IR part of the spectrum to char- 
acterize the gas properties in the nuclear and circum- 
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nuclear regions (Gonzalez-Dclgado et al 1995, here- 
after GD95, and references therein). GD95 conclude 
that the burst of star formation in the nucleus of the 
galaxy should have an age around 4-5 Myr, in view 
of the WR bump detection, and pointed out the pres- 
ence of a previous burst of star formation on the basis 
of the detection of calcium triplet absorption at 8600 
A. Models by Garcia- Vargas et al (1997) indicate that 
a young burst of age between 3.5 and 5 Myr is able 
to explain the emission line spectrum and the WR 
luminosity bump in three circumnuclear HII regions 
of NGC 7714. 

Here, we present HST+GHRS ultraviolet spec- 
troscopy of the nuclear starburst. These data are 
analyzed together with the optical emission line spec- 
trum to obtain information about the ionizing stel- 
lar population. The goal is to perform a consistency 
test between the stellar content derived from multi- 
wavelength observations in the central few hundred 
parsecs in NGC 7714. We will evaluate whether the 
youngest burst can account for the high supernova 
rate (1 yr" 1 ; W81) derived from the X-ray and radio 
continuum emission. A following paper will present 
the results for the NIR K-band spectroscopy and the 
analysis of the multiwavelcngth spectral energy distri- 
bution (SED) (Goldadcr et al 1998; in preparation). 
It will focus on the extincted and low mass stellar con- 
tent of the nuclear starburst in NGC 7714. Calzetti 
(1997) analyzed the properties of the NGC 7714 star- 
burst on a kpc scale. The new GHRS data allow us 
to probe the nuclear starburst on a scale that is an 
order of magnitud smaller. The new data arc useful 
to study if and how the star formation properties in 
NGC 7714 differ between pc and kpc scales. 

The paper is organized as follows. Section 2 
presents the HST observations. Section 3 deals with 
the results and interpretation of the HST ultravio- 
let spectrum. In Section 4 we present photoioniza- 
tion models that fit the emission line spectrum in 
the optical-near infrared spectral range. The mul- 
tiwavelcngth analysis of the NGC 7714 spectrum is 
in Section 5. The summary and the conclusions are 
in Section 6. 

2. HST Observations 

The ultraviolet spectrum of the central nuclear 
starburst in NGC 7714 was obtained in October 1996 
with the Hubble Space Telescope using the GHRS 
and the grating G140L, through the Large Science 



Aperture (LSA, 1.74 x 1.74 arcsec 2 corresponding to 
330 x 330 pc 2 ) . Two different wavelength settings were 
observed, covering 1175 to 1461 A and 1402 to 1689 
A, with a total integration time of 1659 s and 4542 s, 
respectively. After standard pipeline processing, the 
spectra were combined into a single spectrum cover- 
ing 1175 to 1689 A, and rebinned to 0.57 A/pixel, 
which is the nominal dispersion of the grating. The 
spectrum was corrected for redshift assuming a radial 
velocity of 2800 km s _1 (GD95). Figure 1 shows the 
spectrum in the restframe wavelength system of the 
galaxy. 

The spectrum is rich in stellar (CIV A1550, SV 
A1501, CIII A1417, SilV A1400, NV A1240, CIII 
A1175), and interstellar (A1II A1670, Fell A1608, CIV 
A1550, Sill A1527, SilV A1400, CII A1335, OI+SiII 
A1302, Sill A1260) absorption lines. It also shows La 
and Hell A1640 in emission. 

Some of the interstellar lines have weaker satel- 
lites formed in the Milky Way halo (Sill A1527, SilV 
A1400, CII A1335, OI+SiII A1302, and Sill A1260). 
The zero point of the wavelength calibration was 
checked by fitting gaussians to the Galactic absorp- 
tion lines. These lines are blueshifted about 0.3-0.4 A 
with respect to the nominal zero point of -2800 km 
s _1 . This is well within the uncertainties expected 
for an observation through the LSA. The instrumen- 
tal spectral resolution for a point source within the 
LSA is about 0.8 A, but for an extended source it de- 
pends on the size of the ultraviolet source. The widths 
(FWHM) of the gaussians fitted to the Galactic lines 
are 3.3 A or 740 km s _1 at 1400 A. This suggests that 
the object fills a substantial fraction of the aperture. 

The flux measured at 1500 A is 1.9 x 10~ 14 erg 
s _1 cm~ 2 A -1 . This is a factor of 2 lower than the 
IUE flux reported by Kinney et al (1993), indicating 
that the UV bright stellar population significantly ex- 
ceeds the area of the GHRS LSA. An HST WFPC2 
PC image through the filter F606W (A G = 6010.6 A, 
AA= 1497 A) was retrieved from the HST archive. 
The band pass of this filter includes strong emission 
lines such as [OIII]A5007 and Ha, so that the light de- 
tected through this filter is not only continuum emis- 
sion. The integration time was 500 s. The sampling 
is 0.046 arcsec/pixel, corresponding to 8.7 pc/pixel. 
Figure 2a shows a view of the whole galaxy. Most 
of the emission comes from the nucleus with some 
contribution from the two spiral arms. This central 
emission is resolved into a compact knot surrounded 
by several smaller knots. The GHRS spectroscopy 
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suggests that the starburst is extended at ultraviolet 
wavelengths as well. Therefore the spatial morphol- 
ogy in the optical and the ultraviolet is rather similar. 

3. Stellar content in the nuclear starburst in 
NGC7714: fitting the HST UV spectrum 

3.1. Wind absorption lines 

Massive hot stars develop strong stellar winds due 
to radiation pressure in ultraviolet resonance lines 
(Morton 1967; Lucy & Solomon 1970). Typical 
wind velocities in O stars are about 2000 km s _1 
to 3000 km s _1 (Groenewegen, Lamers, & Pauldrach 
1989; Prinja, Barlow, & Howarth 1990). As a re- 
sult, all strong ultraviolet lines in the spectra of O 
stars originate predominantly in the outflow, and have 
blucshiftcd absorptions. The profile shapes reflect the 
stellar mass-loss rates, which are a strong function 
of the stellar luminosity (Castor, Abbott, & Klein 
1975). Since there exists a well-defined stellar mass- 
luminosity relation, the line profiles ultimately de- 
pend on the stellar mass, and — for a stellar popula- 
tion — on the IMF and SF history (Leitherer, Robert, 
& Hcckman 1995). 

The strongest stellar wind resonance features are 
OVI A1034, NV A1240, SilV A1400 and CIV A1550. 
In massive stars, these lines form above the photo- 
sphere in the stellar wind — either as a blueshifted 
absorption in stars with weak winds, or as a P Cygni 
profile if the wind density is sufficiently high. OVI 
and CIV are strong lines in O stars of all luminosity 
classes (Walborn, Bohlin, & Panek 1985; Gonzalez 
Delgado, Leitherer, & Heckman 1997). In contrast, 
SilV is luminosity dependent, and only blue super- 
giant stars produce a strong P Cygni profile (Walborn 
et al 1985). The recombination line Hell A1640 can 
also be formed in very massive O and WR stars with 
very dense winds. Evolutionary synthesis models for 
CIV (Leitherer, Robert, & Heckman 1995) show that 
the line strength increases for a shallower IMF. The 
most massive stars are the dominant contributors to 
the lines, and increasing their numbers relative to the 
less massive stars which produce the continuum - 
either via flattening the IMF slope or increasing the 
upper cut-off mass — leads to stronger lines. SilV 
shows a critical dependence on the age of the burst. 
This spectroscopic method is only efficient in regions 
of recent star formation activity and high stellar den- 
sity where the highest masses can be sampled. The 
technique has been succesfully applied to starburst 



galaxies (Conti, Leitherer, & Vacca 1996; Leitherer 
et al 1996; Gonzalez Delgado et al 1998a) and Seyfert 
galaxies (Heckman et al 1997; Gonzalez Delgado et 
al 1998b). In all these cases, the method has con- 
strained the age and the mass spectrum of the young 
population. 

NGC 7714 has been observed previously by IUE 
(W81; Kinney et al 1993). The spectra show strong 
CIV and SilV in absorption, indicating the presence 
of massive stars in NGC 7714. However, the inter- 
stellar lines are also very strong in NGC 7714 and 
contribute significantly to the strength of the CIV 
and SilV lines. The interstellar contribution can some 
times dominate the profile of the line, as is the case 
of NGC 1705 (Heckman & Leitherer 1997). This is 
why an analysis based only in the equivalent width of 
these lines cannot predict correctly the stellar content 
and the star formation history of the starburst. To 
fully resolve the stellar from the interstellar contribu- 
tion, a resolution better than the IUE resolution of 
R=1000 is required. 

The most conspicuous stellar lines in the GHRS 
spectra of NGC 7714 are the wind lines Hell, CIV, 
SilV and NV. We use line profile synthesis of the ul- 
traviolet wind lines to derive the stellar content and 
to constrain the star formation history of the nuclear 
starburst. We do not fit the entire line profile of the 
CIV and SilV due to the contribution of the interstel- 
lar lines, which are not fully accounted for in the mod- 
els, and are stronger in the starburst spectra than in 
the synthetic models. We select two windows, corre- 
sponding to the blue side of the wind profile (1381.75- 
1390.75 A for SilV and 1528.0-1543.75 A for CIV) and 
to the red side (1403.5-1409.5 A for SilV and 1552.0- 
1560.25 A for CIV), where the profile is dominated 
by the stellar contribution. Then, we compute the 
X 2 parameter between the observations and models. 
The results for SilV are plotted for burst and con- 
tinuous star formation (csf) models as a function of 
the age, upper mass cut-off, and slope of the IMF in 
Figures 3a, b and 3c, d, respectively. The minimum \ 
values for SilV are obtained for instantaneous burst 
models and a Salpeter IMF. The fit for SilV requires 
the presence of O supergiant stars, which appear be- 
tween 3 and 6 Myr after the burst onset. The best 
fit is obtained for a 5 Myr burst. The most massive 
stars present in the starburst have at least a mass 
of 40 MQ. However, we cannot constrain the upper 
mass cut-off because stars more massive than 40 M© 
have already reached the supernova phase at 5 Myr. 
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Thus, the upper mass cut-off is higher than 40 M©. 
Figure 4 shows the profiles of the observed CIV and 
SilV and the 5 Myr burst model (Salpeter IMF and 
M„ p =100 MO). The CIV profile is equally well fitted 
by csf models, but it is not possible to constrain the 
age with these models. 

Note that the emission of the NV P Cygni profile 
is weaker than predicted by the 5 Myr burst model. 
This discrepancy could be due to the Sill A1260 ab- 
sorption line formed in the Milky Way halo that falls 
on the blue edge of the emission part of the profile. 

The determination of the age is also consistent with 
the presence of WR stars in the starburst. WR fea- 
tures have been detected at optical wavelengths (van 
Brcugel et al 1985; GD95), and also in our ultravio- 
let spectrum. The Hell A1640 is a recombination line 
which shows a broad emission profile if it is formed in 
the very dense stellar winds of WR stars and 03-05 
supergiants (Walborn, Bohlin, & Panek 1985). This 
line is prominent in the GHRS spectrum, and the ob- 
served flux is 4x 1CT 14 erg s _1 cm~ 2 . Later we will use 
this flux to estimate the number of WR stars present 
in the nuclear starburst. 

3.2. Photospheric absorption lines 

The spectrum shows several photospheric lines, 
such as SV A1501, CIII A1427, and Silll A1417. These 
lines are prominent in O and early B stars (Walborn 
ct al 1985; Walborn, Parker, & Nichols 1995); they are 
not resonance lines, and therefore they do not form in 
the interstellar medium. Their measured equivalent 
widths are in Table 1. Note that the strength of these 
lines is also consistent with the synthetic model (see 
Figure 5). 

3.3. Extinction Estimates 

To derive the number of WR and O stars in the 
starburst, we need to estimate the extinction. This 
can be done using the ultraviolet continuum flux dis- 
tribution. Lcitherer & Heckman (1995) have derived 
from evolutionary synthesis models that the UV en- 
ergy distribution arising from a starburst has a spec- 
tral index, a, in the range -2.6 to -2.2 (Fa oc A' 3 ), if 
the burst is less than 10 Myr old. This spectral in- 
dex is rather independent of the metallicity and IMF. 
Therefore, any deviation from the predicted spectral 
index can be attributed to reddening. 

First, the spectrum is corrected for the Galac- 
tic extinction, which we derive assuming a ratio of 



N#//E(B-V)= 4.93 x 10 21 cm 2 (Bohlin 1975) and an 
HI column density towards NGC 7714 of 4 x 10 20 
cm~ 2 (Stark et al 1992). The color excess derived 
is E(B — V) = 0.08; this implies an extinction of 
0.65 mag at 1500 A. The slope measured in the cor- 
rected spectrum (log F\ vs. log A) after applying 
the correction for foreground extinction is /3=-2.45, 
very close to the theoretical value of -2.5 which corre- 
sponds to the age deduced from the fit to the absorp- 
tion profiles. This indicates very low additional inter- 
nal extinction associated with the starburst. Using 
the empirical extinction law from Calzetti, Kinney, & 
Storchi-Bcrgmann (1994) to match the slope /3=-2.45 
to the expected slope /3=-2.5 at 5 Myr, we derive an 
internal extinction in the starburst of E(B—V) = 0.03 
(0.26 mag at 1500 A). The spectrum corrected for to- 
tal extinction matches the slope of the synthetic spec- 
trum that fits the wind absorption lines (see Figure 5). 
We note that the total extinction derived from the ul- 
traviolet continuum is lower than the total extinction 
derived from the Balmer decrement, E(B-V)~0.21 
(GD95). This discrepancy between the extinction val- 
ues derived from the two methods has previously been 
found by Fanclli, O'Conncll, & Thuan (1988) and 
Calzetti, Kinney, & Storchi-Bcrgmann (1994). This 
result suggests that probably ionized gas and stars 
are not co-spatial in a starburst (Calzetti 1996). 

The extinction derived above pertains only to the 
ultraviolet light within the GHRS aperture; however, 
the IUE spatial profile indicates that the ultraviolet 
light is extended, at least along P.A.=135°.45. Fig- 
ure 6 shows the IUE ultraviolet light profile that was 
obtained by adding all the data in the dispersion di- 
rection in the wavelength range from 1250 to 1850 
A. The IUE aperture was oriented at P.A= 135°. 45. 
The profile shows a wing on one side that extends be- 
tween 4 and 10 arcsec from the camera pixel with the 
highest emission. This wing could represent the ultra- 
violet emission from the circumnuclear region called 
A by GD95. We have followed the same procedure to 
determine the extinction for the IUE spectrum. After 
correcting for Galactic extinction, the spectral index 
of the corrected IUE spectrum is /3=-l.l. This re- 
sult is in agreement with that found by Heckman et 
al (1998), /3=-1.03. The change of the slope with re- 
spect to the GHRS spectrum indicates a change of 
the extinction and/or of the stellar population with 
the aperture. The change of the extinction with the 
aperture could be due to the swept dust by the stel- 
lar winds and SN explosion. Using the empirical ex- 
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tinction law from Calzetti et al (1994), we derive an 
internal extinction E(B - V) = 0.3 (2.8 mag at 1500 
A). This represents an upper limit to the reddening, 
because an extended population of B and A stars 
could contribute significantly to the ultraviolet con- 
tinuum light in the IUE aperture, flattening the spec- 
tral slope, in which case the intrinsic reddening of the 
starburst would be lower. Thus, the IUE spectrum re- 
veals the presence of extended star formation which 
could be more dusty than the nuclear burst and/or 
more evolved. Alternatively, the IUE aperture could 
pick up a large flux contribution from an older field 
population, that makes the slope of the spectrum flat- 
ter than expected by a hot stellar population. Note 
that the color excess derived depends on the extinc- 
tion law used. If we use the LMC or the SMC law, the 
E(B-V) derived is 0.2 and 0.1, respectively, instead of 
0.3. 

3.4. Intrinsic Luminosities and the WR/O ra- 
tio 

After correcting for the total extinction, E(B-V)= 
0.11, the monochromatic luminosity at 1500 A in the 
GHRS aperture is 10 39 9 erg s _1 A -1 . Assuming that 
the ultraviolet flux is due to a burst of 5 Myr, and 
a Salpeter IMF with upper limit mass cut-off of 100 
MO, the corresponding number of O stars is about 
16600. Note that an upper cut-off mass of 40 M© 
gives the same result. These stars produce a photon 
ionizing luminosity of 7.9 x 10 52 s . We have mea- 
sured the nuclear H/3 flux from the 2D data in GD95 
within two different apertures of sizes that bracket 
the GHRS aperture size, 1.2 x 3.3 arcsec 2 and 1.2 x 2 
arcsec 2 . After correcting for the extinction derived 
from the Balmer decrement, the ionizing photon lu- 
minosity deduced from the R(3 flux within these two 
apertures is 8.2 x 10 52 s" 1 and 6.0 x 10 52 s _1 , re- 
spectively. This is in very good agreement with the 
ionizing photon luminosity deduced from the UV con- 
tinuum. However, some H/3 flux could be lost due to 
seeing and the width of the slit. Therefore a conser- 
vative upper limit to the ionizing photon luminosity 
is 1.4-1.6 times the values derived above. The esti- 
mated mass of the starburst is at least 5 x 10 6 M© 
(integrating the IMF down to a lower mass limit of 

1 M©), and the bolometric luminosity is 5 x 10 9 L©, 
similar to the nuclear starbursts observed in Seyfert 

2 galaxies (Gonzalez Delgado et al 1998b). The ion- 
izing photon luminosity predicted by continuous star 
formation models is 10 53 5 ph s _1 , a factor 4 larger 



than the value deduced from the H/3 flux. 

We can estimate the number of WR stars from the 
luminosity of the Hell A 1640 line measured within 
the GHRS aperture. After dereddening, and assum- 
ing the calibration given by Conti (1996; 1 WN cor- 
responds to 6.3 x 10 36 erg s _1 ), the line flux gives a 
luminosity of 10 40 1 erg s _1 , corresponding to about 
2000 WR stars. This number of WR stars gives a 
WR/O ratio of 0.12. The equivalent width of the Hell 
A1640 is 2.7 A. These two numbers are compatible 
with the values predicted by an instantaneous burst 
(see Figures 2 and 10 in Schaerer & Vacca 1998). Our 
conclusion is that a burst 5 Myr old with a mass of 
5xl0 6 M© is able to explain the continuum and stel- 
lar UV absorption features in the central 1.74x1.74 
arcsec of NGC 7714, with a Salpeter IMF and upper 
mass limit higher than 40 MQ. 

The monochromatic luminosity at 1500 A in the 
IUE aperture is 10 41 25 erg s _1 A -1 , after correct- 
ing for a total reddening E(B-V)=0.38. This implies 
a bolometric luminosity of 1.1 x 10 11 LQ, an ioniz- 
ing photon luminosity of 2xl0 54 ph s" 1 and a mass 
of 1.2 x 10 s M© (assuming an age of 5 Myr), a fac- 
tor 20 larger than the values derived from the GHRS 
aperture. The IUE bolometric luminosity is similar 
to the total IR luminosity, 5 x 10 10 LQ, measured by 
IRAS and computed with the four wavelength bands 
as defined by Sanders & Mirabel (1997). Thus, this 
extended star-formation region could contribute sig- 
nificantly to the overall galaxy bolometric luminosity. 
However, these results have been derived under the 
assumption that the change in the slope in the IUE 
spectrum is due to internal extinction. If an under- 
lying older population contributes significantly to the 
IUE aperture by flattening the spectrum, the derived 
color excess, and the mass, bolometric luminosity and 
the ionizing photon luminosity, represent an upper 
limit to the true values. The other extreme case is 
to assume that no internal extinction affects the IUE 
flux. In this case, the monochromatic luminosity at 
1500 A is 10 40 1 erg s _1 A -1 . This implies an ioniz- 
ing photon luminosity of 1.2x10 ph s _1 , which is 
a factor 2.5 lower than the rate of ionizing photons 
derived from the Ha image adding the contribution 
of the nucleus and the region A (note that these two 
components fall within the IUE aperture, see Figure 1 
in GD95). Thus, the change of slope is probably due 
to both, a change of reddening and the contribution 
of an underlying population. 

Evidence that an underlying population contributes 
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to the IUE aperture comes from the CIV profile. This 
line is significantly more diluted in the IUE than 
in the GHRS spectrum (Figure 7). Thus, the IUE 
spectrum suggests that an underlying population con- 
tributes to its flux, and this makes the determination 
of the starburst properties from the IUE flux rather 
uncertain. 

3.5. Interstellar lines and La emission 

The ultraviolet spectrum shows two systems of in- 
terstellar lines, one formed in the Milky Way halo 
and the other in NGC 7714 (see Figure 1). The 
most prominent lines in NGC 7714 are A1II A1670, 
Fell A1608, Sill A1527, CII A1335, OI+SiII A1302, 
Sill A1260 and Sill A1190. The equivalent widths of 
these lines are shown in Table 1. Their values are 
uniformly distributed between 2 and 3 A, indicating 
that the lines are saturated. In the saturated part of 
the curve of growth the line equivalent width depends 
on the velocity dispersion rather than on the column 
density, so it can be used to constrain the kinemat- 
ics of the interstellar medium (Hcckman & Leitherer 
1997; Gonzalez Delgado et al 1998a). The equivalent 
width of the Sill A 1260 implies a velocity dispersion 
of 160 km s _1 , which is consistent with the upper 
limit derived from the FWHM of these lines. Note 
that the lines are not resolved since our instrumental 
resolution is 3.3 A, corresponding to a velocity disper- 
sion of about 300 km s _1 . The observed equivalent 
width could be the result of large scale motions of 
the interstellar gas. This has also been observed in 
other starburst galaxies (Heckman & Leitherer 1997; 
Gonzalez Delgado et al 1998a). However, in contrast 
to other starbursts, the ultraviolet absorption lines in 
NGC 7714 are not blueshifted with respect to their 
systemic velocity. 

La is observed in emission. The measured flux is 
2.9 xl0~ 14 erg s _1 cm~ 2 and the equivalent width 
is 1.5 A. The line shows a pronounced drop in the 
blue wing, due to interstellar La absorption in our 
Galaxy and in NGC 7714 itself. This absorption 
shifts the peak of the emission line to the red with re- 
spect to the systemic velocity. This has also been ob- 
served in other starburst galaxies (Lequeux et al 1995; 
Gonzalez Delgado et al 1998a; Kunth et al 1998). 
To correct the La emission from the absorption, we 
have performed a multiple component fit to the ob- 
served absorption features (Sill A1190,1193, Galactic 
La and La in NGC 7714; we assume that La in NGC 
7714 is at the rest wavelength) using theoretical Voigt 



profiles produced with the XVoigt software package 
(Alar & Bailey 1995). After dividing the observed 
spectrum by the fitted Voigt profiles, the La flux is 
2.9 x 10~ 13 erg s" 1 cur 2 and the La/U(3 ratio 2.9 
(4.1 using the H/3 flux measured in the central 1.2x2 
arcsec 2 ). Correcting by the reddening derived from 
the Balmer decrement and using the LMC extinction 
law, the ratio is 23 (33); close to the value predicted 
from the recombination theory (33, Ferland & Oster- 
brock 1985). We can conclude that attenuation by 
dust due to multiple resonant scattering by hydrogen 
atoms does not play an important role in NGC 7714 
because the La/H/3 ratio corrected for extinction is 
close to the recombination value. 

4. Stellar content in the nuclear starburst: 
fitting the optical-near infrared emission 
line spectrum 

The emission line spectrum of a starburst depends 
on the radiation field from the ionizing stellar clus- 
ter, and the density and chemical composition of the 
gas. In this section, we predict the emission line spec- 
trum using a photoionization code that takes as in- 
put the spectral energy distribution generated by a 
stellar evolutionary synthesis code. This technique is 
applied as a second independent method to constrain 
the star formation history and the age of the nuclear 
starburst in NGC 7714. It has been previously used to 
study the stellar content in starbursts and giant HII 
regions (Garcia- Vargas & Diaz 1994; Garcia- Vargas, 
Brcssan, & Diaz 1995a, 1995b; Stasihska & Leitherer 
1996; Garcia- Vargas et al 1997). The goal is to per- 
form a test of consistency between the young stellar 
content derived from the ultraviolet absorption spec- 
trum and from the optical emission lines. First, we 
describe the basis of the stellar evolutionary synthesis 
and photoionization models. 

4.1. Stellar evolutionary synthesis models 

The spectral energy distribution was generated by 
the evolutionary synthesis code developed by Lei- 
therer and collaborators (Leitherer, Robert, & Dris- 
sen 1992). The code has been updated recently. A 
description of the new version will be presented else- 
where (Leitherer et al 1998, in preparation). Two im- 
portant changes have been made with respect to the 
previous version of the code. The most recent stellar 
evolutionary models of the Geneva group (Schaller et 
al 1992; Maeder 1994), and the stellar atmospheres 
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grid compiled by Lejeune et al (1996) are used, in- 
stead of Maeder's (1990) and Kurucz's models (1992). 
As in the previous version of the code, we use the 
expanding spherical extended non-LTE models pub- 
lished by Schmutz et al (1992) for stars with very 
strong winds. 

The spectral energy distribution was generated us- 
ing the 0.4 Z© metallicity tracks (note that the chem- 
ical abundance derived from the emission lines in 
NGC 7714 is close to half solar, GD95), and assum- 
ing two different star formation scenarios (instanta- 
neous burst and continuous star formation). In both 
cases, a Salpeter IMF with an upper and lower mass 
limit cut-off of 80 M© and 1 MQ, respectively, was 
assumed. 

4.2. Photoionization models 

We take the spectral energy distribution as input 
to the photoionization code CLOUDY (version 90.04, 
Ferland 1997). CLOUDY resolves the ionization- 
recombination and heating-cooling balances, and pre- 
dicts the ionization structure of the nebula, the tem- 
perature, and the intensities of the emission lines. We 
assume that the nebular gas is ionization bounded 
and spherically distributed around the ionizing clus- 
ter with a constant density. We assume an inner ra- 
dius of 3 pc, but the outer radius is determined by 
the ionization front. The chemical composition of the 
gas is scaled to the oxygen abundance (close to half 
solar), except for He, S, Ne and N. The abundances 
and the electron density were determined from the 
optical emission lines as given by GD95. They are 
shown in Table 2. The adopted values of the solar 
abundance are as in Stasihska (1990), except the sul- 
fur abundance, for which we take the S/O ratio given 
by Grevesse & Anders (1989). The ionizing photon 
luminosity is fixed to log(Q)=52.9 ph s^ 1 , as derived 
above independently from the Balmer recombination 
lines and from the ultraviolet continuum luminosity. 
Models are computed taking the filling factor as a free 
parameter, with values of 10 _1 , 10~ 2 , 10~ 3 , 10~ 4 , and 
10 -5 . The change in the filling factor is equivalent 
to changing the ionization parameter U, defined as 
Q/(47rR s N e c); where Q is the ionizing photon lumi- 
nosity, N e the electron density, c the speed of light 
and R s the Stromgren radius. The average U is pro- 
portional to (<j) 2 N e Q) 1 / 3 , where 4> is the filling factor. 

The emission line spectrum depends on the ion- 
izing radiation field — i.e. on the evolution of the 
cluster — and on the ionization parameter for a fixed 



geometry, metallicity, density and the total number 
of ionizing photons. First, we determine the ioniza- 
tion parameter using the ratio [SIIJA6716+6731/H/3. 
This ratio is a good calibrator of U for continuous star 
formation and for burst models (Figure 8). The ob- 
served ratio [SII]A6716+6731/H/3=0.58 indicates an 
ionization parameter of -2.9, corresponding to a fill- 
ing factor of 0.001. Then, we use the other emission 
lines to constrain the star formation scenario (burst 
or continuous star formation) and the age of the star- 
burst. Figure 9 shows the line ratios [OIIL1A5007/H/3, 
[SIII]A9069/H/3, [OII]A3727/H/3, [NII]A6584/H/3, [OI]A6300/H/3 
and HcIA4471/H/3 as a function of age for continous 
star formation and burst models with filling factor 
0.001. Burst models predict line ratios that arc in 
better agreement with the observed ratios than con- 
tinuous star formation models. Table 3 shows the 
strength of the emission lines observed and the val- 
ues predicted by the 4.5 Myr burst. The 4.5 Myr 
burst model fits well the strongest emission lines, but 
it does not reproduce the [OIII]A4363 and [OI]A6300 
lines; the values observed are a factor of two larger 
than predicted. GD95 did three different estimations 
of the electron temperature of the nebula using the 
ratio of the auroral to the nebular collisionally ex- 
cited lines [SIII], [SII] and [Nil] lines. The three ra- 
tios give a similar electron temperature (about 7500 
K), and with this electron temperature the metallic- 
ity was derived. However, the observed [OIIIJA4363/ 
[OIII]A5007 ratio is larger than the expected ratio for 
a nebula with a electron temperature of 7500 K. This 
is probably a consequence of a larger contribution of 
the [OIIIJA4363 line to the temperature ratio. In a 
nebula the presence of shocks can significantly af- 
fect the lines [OIII]A4363 and [OI]A6300 (Peimbert, 
Sarmiento, & Fierro 1991). Shocks are naturally ex- 
pected to be present in the ionized gas surrounding 
a starburst with age of 5 Myr as a result of the stel- 
lar wind and the explosion of the first supernovae in 
the stellar cluster. At this age, the deposition rate of 
mechanical energy is about 10% of the ionizing lumi- 
nosity. Thus, shock heating could probably account 
for the excess emission in [OIIIJA4363 and [OIJA6300 
with respect to the photoionization values. 

This photoionization model predicts a Stromgren 
radius of 185 pc. If the nebula is radiation bounded, 
the radius is a function of the ionizing photon flux, 
the electron density and the filling factor. Then for 
fixed Q, Ne and the radius has to be equal to the 
size of the nebula. We know that the ionizing cluster 
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is not a point source because the ultraviolet interstel- 
lar absorption lines in the GHRS spectrum are not 
resolved. It is at least as extended as the size of the 
GHRS aperture (330x330 pc). This size is similar 
to the Stromgren radius predicted. Also, from the 
WFPC2 image (see Figure 2), we estimate that the 
size of the nebula is about 210 pc. Thus, the predicted 
size of the Stromgren radius is compatible with the es- 
timated size of the HII nebula. This result is also con- 
sistent with the HII region being radiation bounded. 
Further evidence comes from the [OIJA6300 line. This 
low ionization line forms in the partially ionized zone, 
which does not exist in a matter bounded nebula. On 
the other hand, the ionizing photon flux derived from 
the ultraviolet continuum luminosity is very similar to 
that derived from the H/3 luminosity, indicating again 
that the radiation bounded hypothesis is correct. 

Thus, the conclusion is that emission line ratios can 
discriminate very well between the burst and contin- 
uous star formation (csf) scenario. The csf models 
can be excluded because they predict emission line 
ratios which are higher than the values observed (ex- 
cept for the [Nil] A6584/H/3). The best fit is found 
for a burst 4.5 Myr old. Is it reasonable that all these 
stars are formed in an instantaneous burst? The ul- 
traviolet data indicate that the starburst is extended 
at least 1.74 arcsec (330 pc). The HST optical image 
shows that it may not be more extended than 400 
pc (in diameter). Taking 400 pc as the diameter of 
the nuclear starburst and the velocity dispersion for 
this galaxy of 180 km s _1 (GD95), we derive that the 
crossing time within the starburst is 2.1 Myr. This 
time is significantly shorter than the evolutionary life- 
time of a single generation of massive stars. Thus, it 
is dynamically possible that the nuclear starburst in 
NGC 7714 formed in an instantaneous burst of star- 
formation because the crossing-time is shorter than 
the duration of the starburst. 

5. Discussion 

5.1. X-ray and radio emission 

Can the derived stellar population explain the ob- 
served X-ray luminosity and radio continuum emis- 
sion? X-ray emission in normal galaxies can have 
several origins. As discussed by Martin & Kennicutt 
(1995) and Stevens & Strickland (1998), production 
of X-rays in starburst galaxies by AGNs, static galac- 
tic coronae, single stars, or X-ray binaries is rather 
unlikely to be the dominant mechanism. The most 



plausible interpretation is due to powering by super- 
novae, with some contribution from stellar winds. 

The models of Leitherer & Heckman (1995) pre- 
dict a type II supernova rate of about 0.007 yr~ 4 
for a bolomctric luminosity of 5 x 10 9 L0. If we 
assume that the kinetic energy released per super- 
nova is 10 51 erg, the mechanical luminosity released 
by stellar winds and supernovae is 2 x 10 41 erg s _1 . If 
thermalization occurs, this luminosity is in principle 
available to create a hot, over-pressured bubble. In 
practice, the sum of the thermal and kinetic energy of 
the bubble will be lower due to radiation losses. Hy- 
drodynamical simulations by Thornton et al. (1998) 
suggest radiation losses on the order of 90%. The ob- 
served X-ray luminosity in the center of NGC 7714 is 
6 x 10 40 erg s^ 1 (W81). If interpreted as due to ther- 
mal emission from hot, shocked material, it could be 
accounted for by supernovae and winds, even allowing 
for significant radiation losses. 

The observed 20 cm radio flux is 20 mJy (W81). 
In the standard supernova driven outflow model (e.g., 
Chevalier 1982; 1991) radio emission is assumed to 
be synchrotron radiation from electrons accelerated 
by shocks in the interaction zone. Models for the 
non-thermal component of a hot bubble are substan- 
tially more uncertain than the thermal (X-ray) pre- 
diction. Empirical data for individual local super- 
nova remnants indicate non-thermal luminosities of 
about 1% the thermal value (Chevalier 1982). As- 
suming the 20 cm flux is purely non-thermal, and 
adopting a bandwidth of Af = 10 GHz, the non- 
thermal component in NGC 7714 has a luminosity 
of 3.7 x 10 38 erg s -1 . Remarkably, this is about 1% 
of the X-ray luminosity. Therefore we conclude that 
both the thermal and the non-thermal component of 
the hot, central ISM in NGC 7714 are consistent with 
being powered by winds and supernovae. 

While simple energy considerations are instructive, 
more detailed hydrodynamical calculations for the be- 
havior of a hot supcrbubble power by a stellar pop- 
ulation arc needed. Although still deficient in some 
aspects (e.g. with respect to the time dependence of 
the energy input), the relations discussed by Mar- 
tin & Kennicutt (1995) and Stevens & Strickland 
(1998) provide further insight. These authors used 
the standard superbubble model to solve the energy 
and momentum equations for an adiabatic bubble. 
According to these models, the X-ray luminosity can 
be expressed in terms of the starburst luminosity and 
age for a given IMF. A 5 Myr old starburst with a 
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standard Salpeter IMF is predicted to have a ratio 
of the X-ray over bolometric luminosity of ~ 10 -3 
(see Fig. 16 of Stevens & Strickland). Applied to 
NGC 7714, these models predict an X-ray luminos- 
ity of 2 x 10 40 erg s _1 , within a factor of 3 of the 
observed value. The fact that the prediction is lower 
than the observation might suggest a clumpy ISM. An 
ambient medium containing clouds or filaments can 
generate extra mass, and therefore extra X-ray lumi- 
nosity for fixed mechanical luminosity input (Martin 
& Kennicutt). However, the theoretical and obser- 
vational uncertainties are large enough that the data 
are consistent with a homogeneous medium as well. 

The stellar population, as derived independently 
from nebular emission lines and stellar absorption 
lines, is capable of generating the mechanical energy 
input to account for the observed thermal X-ray and 
non-thermal radio emission as well. The associated 
supernova rate is 0.007 yr _1 . This rate is consistent 
with the upper limit imposed by the null result in the 
supernova search of Richmond, Filippenko, & Galisky 
(1998). No supernovae were detected in NGC 7714 
during a several year long monitoring period. The 
non-detection implies an upper limit to the supernova 
rate of < 0.1 yr -1 . Few if any supernovae should es- 
cape detection due to dust extinction since the extinc- 
tion is known to be low. We can account for most of 
the luminosity at 1500 A, which traces precisely the 
progenitors of type II supernovae (M w 20 MQ). 

The derived supernova rate is two orders of mag- 
nitude lower than that found by W81. Their esti- 
mate was based on using X-ray and radio luminosi- 
ties and lifetimes of local supernova remnants as tem- 
plates and scaling them to the luminosities observed 
in NGC 7714. In retrospect this approach has turned 
out to be too simple. Supernovae and supernova rem- 
nants in dense media, like in starbursts, can have 
rather different properties. This may be particularly 
relevant to the very young starburst in NGC 7714. 
Few red supergiants have formed at an age of 5 Myr, 
and most of the supernova progenitors are predicted 
to be WR stars. WR stars have hydrogen-free en- 
velopes and have stellar-wind properties quite differ- 
ent from red supergiants. On the other hand, W81 use 
the truncated starburst models of Riekc ct al (1980) 
to estimate the mass and rate of supernovae. In these 
models the IMF is truncated with the high mass < 
30 MQ. Therefore, they have to boost the star for- 
mation rate in oder to get the ionizing photons, and 
therefore the supernova rate became high. However, 



our data are compatible with a normal IMF. 

In the last few years, young compact supernovae 
have been studied in detail at radio wavelengths 
(Weiler et al, 1986, 1990, 1991, 1992; Van Dyk et 
al 1992, 1993a,b; Ryder et al 1993; Yin 1994). These 
supernova, at the peak of their emission, can have 
a radio luminosity which is 100 - 1000 times larger 
than Cas A. The presence of these supernovae and 
their remnants can explain the radio luminosity de- 
tected in highly luminous IRAS galaxies (Colina & 
Perez-Olea, 1992; Perez-Olea & Colina, 1995). Evi- 
dence that these compact supernova can account for 
the radio emission comes from 18 cm VLBI contin- 
uum imaging observations of Arp 220. The images 
have revealed a dozen compact radio sources, with 
fluxes that can be explained by compact young radio 
supernova (Smith et al 1998). 

Assuming that the 6 cm radio emission observed 
by W81 is non-thermal emission emitted by the cen- 
tral 1.74x1.74 arcsec, and using the models of Colina 
& Perez-Olea (1995), we derive a supernova rate of 
0.07 yr _1 , which is in agreement with our empirical 
upper limit determination, but still one order of mag- 
nitude higher than the rate predicted by the model 
which accounts for the UV luminosity. To predict a 
supernova rate as high as 0.07 yr _1 , we need a cen- 
tral starburst with a UV luminosity higher than 10 39 - 9 
erg s^ 1 A -1 . Is the estimated luminosity wrong be- 
cause of a wrong estimation of the intrinsic redden- 
ing derived from the UV continuum slope? If the 
internal reddening were E(B — V) = 0.25, then the 
UV luminosity of the starburst in the GHRS aper- 
ture would be 10 40 - 7 erg s _1 A -1 , the supernova rate 
0.05 yr _1 , the thermal emission 4 mJy, and the non- 
thermal emission 6 mJy. However, by correcting the 
GHRS spectrum by E(B - V) = 0.25, the UV spec- 
trum appears clearly over-corrected when compared 
with the synthetic model. From exploring the fit of 
the UV continuum slope, we estimate that the intrin- 
sic reddening cannot be higher than 0.1. Correcting 
for this extinction, the luminosity at 1500 A is 10 40 ' 2 
erg s^ 1 A -1 , the mass of the starburst 1.1 xlO 7 M© 
and the supernova rate 0.02 yr _1 . In this case, we 
can account for just one third of the radio emission. 
Compact supernova remnants are also very strong X- 
ray sources with luminosities in the range 10 38 -10 
erg s" 1 (e.g. SN 1978K, Ryder et al 1992; SN 1986J, 
Bregman & Pildis 1992; SN 1993J, Schlegel 1994). 
SN 1988Z was detected from ROSAT HRI observa- 
tions with a luminosity of 10 41 erg s _1 (Fabian & 
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Terlevich 1996). Thus, these compact supernova rem- 
nants arc brighter by several orders of magnitude than 
the Galactic remnants (10 35 -10 37 erg s _1 range), and 
they can contribute significantly to the X-ray flux. 

Probably even more important, a starburst nucleus 
is not just a scaled-up version of an individual super- 
nova remnant. The superposition of numerous super- 
nova explosions and the effects of stellar winds cre- 
ate a hot, over-pressured bubble which expands and 
can experience a break-out from the galaxy (Heckman 
1995). The relation between the X-ray luminosity of 
the bubble and the properties of the stellar popula- 
tion, including the supernova rate, can be quite dif- 
ferent from what is expected on the basis of individ- 
ual supernova remnants (although the physics is re- 
lated) . The success of the simple superbubble model 
discussed above is encouraging. We conclude that the 
previously derived supernova rate of ~1 yr _1 is not 
supported by more recent starburst models and that 
the current supernova rate in the nucleus of NGC 7714 
is close to about 1 per century. 

Incidentally, this rate is about the same as the disk- 
integrated rate observed in normal spirals, including 
the Milky Way. NGC 7714 does not have a signif- 
icantly higher total supernova rate than observed in 
spirals with active star formation. What distinguishes 
NGC 7714 (and other nuclear starbursts), is the con- 
centration of the supernova events in the nucleus. The 
nuclear supernova rate per unit surface is about three 
orders of magnitude higher than that of normal galax- 
ies. 

5.2. Evidence of an additional stellar popula- 
tion 

Is there evidence for an additional stellar compo- 
nent younger than 30-40 Myr in the nucleus of NGC 
7714? Bursts younger than 30-40 Myr will contribute 
to the non-thermal emission (stars more massive than 
8 M© explode as supernovae, and their life-time is 
less than 30-40 Myr). When older than 10 Myr, they 
will not contribute to the strength of the emission 
lines, but will affect the equivalent widths of the re- 
combination emission lines. For a burst 4.5 Myr old, 
the equivalent width of H/3 is 65 A, however, the ob- 
served value is 30 A. Thus, the equivalent width of 
H/3 must be diluted. But we do not know whether 
it is diluted by a burst younger than 30-40 Myr, or 
by the bulge population. In the near-infrared, the 
spectrum shows the Call triplet in absorption. The 
strength of these lines is dependent on both metallic- 



ity and gravity (Diaz, Terlevich, & Terlevich 1989). 
At solar metallicity, these lines are very strong in star- 
bursts if the near-infrared light is dominated by red 
supergiant stars. Population synthesis models of the 
two strongest absorption lines (Call A8542 and Call 
A8562) for half solar metallicity predict an equivalent 
width for the two lines of 4-5 A for burst younger 
than 40 Myr (Garcia- Vargas et al 1998). The ob- 
served value is 4.5 A (GD95), in agreement with the 
prediction of the models. However, the same mod- 
els predict an equivalent width somewhat larger than 
5 A for a a burst older than 100 Myr. Thus, even 
though the strength of the Call triplet is compatible 
with a burst younger than 40 Myr, we cannot ex- 
clude that the lines could be dominated by red giant 
instead of red supergiant stars. Thus, the equivalent 
width of H/3 and the strenght of the Call triplet show 
evidence that an underlying population contributes 
to the central 1.7x1.7 arcsec spectrum of NGC 7714. 
Further evidence comes from the NIR spectrum. We 
will present a full analysis of the UV-optical-NIR con- 
tinuum in paper II (Goldader et al 1998, in prepara- 
tion). 

6. Summary and conclusions 

The main goal of this work is to constrain the most 
recent star formation history in the central 300 pc of 
NGC 7714. To archieve this goal, we observed the nu- 
cleus of NGC 7714 from the UV to the near infrared. 
Consistent with previous observations, our data show 
that NGC 7714 exhibits the spectral dichotomy typ- 
ically observed in starburst galaxies. The ultraviolet 
spectrum is dominated by absorption features formed 
in the intestellar medium and stellar winds, while 
the optical spectrum is dominated by emission lines 
formed in the photoionized interstellar medium. We 
constrain the most recent star formation history in the 
nucleus of NGC 7714 by applying two different tech- 
niques, one based on the ultraviolet spectrum, and 
the other on the optical to near infrared emission line 
spectrum. 

Evolutionary synthesis line profiles of the CIV and 
SilV wind lines are used to constrain the duration of 
the starburst, the age of the stellar population, and 
the IMF parameters (slope and upper limit). We find 
that the profiles of the wind lines are compatible with 
a single burst, an age of 5 Myr and a Salpeter IMF 
with upper mass cut-off higher than 40 M©. Argu- 
ments of causality suggest that is dynamically pos- 
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sible that the nuclear starburst formed an instanta- 
neous burst since the dynamical crossing time within 
the starburst is lower than the evolutionary life time 
of massive stars. Evolutionary synthesis models cou- 
pled to the photoionization code CLOUDY are also 
used to generate the emission line spectrum. Con- 
tinuum star formation models are not consistent with 
the strength of the emission line ratios. The best fit to 
the emission lines is found for a 4.5 Myr burst model. 
From the ultraviolet continuum luminosity we derive 
a mass of the starburst of 5-11 xlO 6 M© (note that 
this probably a lower limit to the mass of the starburst 
because it has been derived from the UV continuum 
luminosity) and a bolometric luminosity of 5-10 xlO 9 
L0, assuming stars down to 1 M© are formed. Very 
little extinction is associated with this nuclear star- 
burst. WR features at ultraviolet (Hell A1640) and 
optical (Hell A4686) wavelengths are detected. Their 
luminosities imply a WR/O ratio of ~ 0.1, which is 
also compatible with the predicted values from the 
burst models. 

This model predicts a supernova rate of 0.007-0.02 
yr -1 , consistent with the rate derived from the radio 
and the X-ray luminosities, 0.05 yr _1 (if the emission 
accounts for the contribution of radio and X ray com- 
pact supernova). This rate is in agreement with an 
upper limit of 0.1 yr _1 based on the non-detection 
of supernova events in the past several years. This 
supernova rate is two order of magnitude lower than 
previously suggested. The supernova rate in the NGC 
7714 nucleus is similar to disk-integrated rates in nor- 
mal spiral galaxies. 
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Table 1 

Equivalent width (A) of the absorption and emission features 3. 



Line 


type 


Ew (A) 


NV A 1240 


W 


2.8 


NV A 1240 


W 


-1.2 


Sill A1260 


IS 


1.7 


Ol+Sill A1300 


IS 


3.1 


CII A1335 


IS 


2.8 


SilV A1400 


W+IS 


6.4 


SilV A1400 


w 


-0.2 


Silll A1417 


Ph 


0.6 


CIII A1428 


Ph 


0.5 


SV A1501 


Ph 


0.8 


Sill A1526 


IS 


2.1 


CIV A1550 


W+IS 


9.4 


CIV A1550 


w 


-0.9 


Fell A1608 


IS 


1.8 


Hell A1640 


w 


-1.6 



a IS means Interstellar, W wind, and 
Ph Photosphcric line. Negative values 
are lines in emission. 



Table 2 

Input parameters for CLOUDY a 



parameter 


value 


log N e 


2.75 


12+log O/H 


8.53 


log He/H 


-1.00 


log S/H 


-5.05 


log Ne/H 


-4.30 


log N/H 


-4.15 



a The electron density 
and metallicity are from 
GD95, and they were de- 
rived from the optical 
emission lines. 
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Table 3 

Line ratio of the emission lines with respect to H/3 a 



Line 


observed 


model 


[Oil] A3727 


1.97 


2.13 


[OIII] A4363 


0.01 


0.005 


Hel A4471 


0.037 


0.04 


[OIII] A5007 


1.59 


1.83 


Hel A5876 


0.12 


0.12 


[SO] A6312 


0.007 


0.01 


[Ol] A6300 


0.043 


0.022 


[Nil] A6584 


1.22 


1.37 


[SII] A6716 


0.29 


0.32 


[SII] A6732 


0.28 


0.32 


[SIII] A9069 


0.27 


0.32 



a The observed line ratios are from 
GD95 corrected for reddening with 
C(H/3)=0.3. The line ratios are pre- 
dicted using the spectral energy dis- 
tribution of the 4.5 Myr burst model, 
assuming a Salpctcr IMF and M up = 
80 M0, as the ionizing radiation field 
into CLOUDY. 
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Fig. 1— The 1.74x1.74 arcscc 2 GHRS spectrum of 
the nucleus of NGC 7714 (in the rest frame of NGC 
7714 wavelength). The most important wind, photo- 
spheric and interstellar absorption lines are labelled. 
The interstellar absorption lines formed in the halo of 
our Galaxy are also labeled as Gal. 

Fig. 2. - The WFPC2 F606W image. The orienta- 
tion is North up and East to the left. The scale is 
0.046 arcsec/pix. (a) The image shows the morpho- 
logical characteristics of almost the whole galaxy, (b) 
The central 2.5x2.5 arcsec 2 . The nucleus is clearly 
extended, and resolved into severals knots. 

Fig. 3. — x 2 parameter of the fits to the profile of 
SilV for burst models (a,b), and for the continuous 
star formation models (c,d) with different M up (a,c) 
and IMF slope (b,d). The size of the bubble is pro- 
portional to x 2 , with a larger size indicating a worse 
fit. The scale in the small boxes indicates the value 
of x 2 - 

Fig. 4. — The observed profiles (dashed line) of CIV 
(a) and SilV (b) and the 5 Myr burst model (thick 
line), assuming a Salpeter IMF and M up = 100 M©. 



Fig. 5.— The GHRS spectrum dereddened by E(B- 
V)=0.08 using the MW extinction law, and by E(B- 
V)=0.03 using the Calzetti et al (1994) extinction 
law, and the synthetic 5 Myr burst model (in rela- 
tive units). The IMF slope is Salpeter and M up = 
100 M0. 

Fig. 6. — The IUE intensity profile perpendicular to 
the dispersion direction integrating from 1250 to 1850 
A. The zero in the X-axis is the IUE camera pixel 
with highest ultraviolet emission. The bump between 
4 and 9 arcsec could be the emission associated with 
the circumnuclear region called A in GD95. 

Fig. 7. — The GHRS (thick line) spectrum rebinned 
to the resolution of the IUE spectrum (thin line). The 
GHRS flux at 1500 A is also scaled to match the IUE 
flux. The spectra are corrected only by Galactic ex- 
tinction. 

Fig. 8. — Predicted 

emission line ratio [SII]A6716+6732/H/3 as a function 
of the average ionization parameter generated with 
CLOUDY, taking as input the spectral energy distri- 
bution from the evolution of a burst (open symbols) 
and continuous star forming regimes (full symbols). 
The observed value is indicated by a horizontal line. 
See text for details. 

Fig. 9. — Predicted emission line ratios as a function 
of age. The filling factor is 0.001. The observed values 
are indicated by a horizontal line. These ratios are 
compatible with a burst 4.5 Myr old, and they exclude 
the continuous star formation models. 
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